Pulsar Magnetospheric Emission Mapping: 
Images and Implications of Polar-Cap Weather 
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ABSTRACT 

The beautiful sequences of "drifting" subpulses observed in some radio pulsars have been 
regarded as among the most salient and potentially instructive characteristics of their emission, 
not least because they have appeared to represent a system of subbeams in motion within the 
emission zone of the star. Numerous studies of these "drift" sequences have been published, 
and a model of their generation and motion articulated long ago by Ruderman & Sutherland 
(1975); but efforts thus far have failed to establish an illuminating connection between the drift 
phemomenon and the actual sites of radio emission. Through a detailed analysis of a nearly 
coherent sequence of "drifting" pulses from pulsar B0943+10, we have in fact identified a system 
of subbeams circulating around the magnetic axis of the star. A mapping technique, involving 
a "cartographic" transform and its inverse, permits us to study the character of the polar-cap 
emission "map" and then to confirm that it, in turn, represents the observed pulse sequence. On 
this basis, we have been able to trace the physical origin of the "drifting-subpulse" emission to 
a stably rotating and remarkably organized configuration of emission columns, in turn traceable 
possibly to the magnetic polar-cap "gap" region envisioned by some theories. 



Subject headings: 
nonthermal 



MHD — plasmas — pulsars: general, individual (B0943-I-10) — radiation mechanism: 



Since their discovery by Drake & Craft (1968), 
"drifting" subpulses have fascinated both ob- 
servers and theorists and have often been regarded 
as one of the "keys" to understanding pulsar ra- 
dio emission. Without any adequate model, early 
studies of prominent "drifters" {e.g., Taylor & 
Huguenin (1971); Backer (1973)) focussed on de- 
lineating the drift phenomenon. Then, Ruderman 
& Sutherland's (1975) model not only articulated 
a qualitative picture of "drifting" subpulses as a 
system of "carousel beams" circulating around the 
magnetic axis within the (magnetic) polar emis- 
sion region, but they also attributed the rotation 
physically to ExB drift and estimated how long 



these circulation times would be. 

While fascinating and theoretically compelling, 
the phenomenon has proved difficult to interpret 
physically. No such sequence has appeared precise 
enough to determine the structure and dynamics 
of the beams, and thus to see if it follows any pat- 
tern which might or might not confirm the expec- 
tations of the aformcntioned model. Any sightline 
to a star can cross the periphery of its polar emis- 
sion cone at only one or two points, and it was 
difficult to fix accurately just where these points 
fell in relation to the rotational and magnetic axes. 
Nor was it clear to what extent the "drift" rates 
were aliased by the rotation frequency of the star. 
Within the Ruderman & Sutherland model, the 
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total number of subbcams must be ascertained, or 
no calculation of the circulation time is possible — 
and thus no quantitative connection can be made 
with theories of the polar-cap "gap" region. It is 
a telling statement on this unfortunate state of af- 
fairs that the most extensive studies of a pulsar 
with a highly regular "drifting" -subpulse pattern, 
B0809-t-74, were published 25-30 years ago! [Tay- 
lor et al. (1971); Manchester et al. (1975)]. 

Pulsar B0943-I-10 is another pulsar in the same 
class and is remarkable in other ways as well. One 
of a handful of pulsars discovered below 300 MHz 
(Vitkevich, et al. 1969), it weakens markedly at 
higher frequencies and has never been detected 
above 600 MHz. Its profile evolution (and "drift- 
ing" subpulses) identify it as a member of the 
conal single {3^) class (Rankin 1993b), and we 
now know that its steep spectrum is due to a 
sightline traverse so peripheral that it ultimately 
misses the ever narrower high frequency cone. 
The pulsar's long sequences of "drifting" subpulses 
have prompted several different studies [Taylor & 
Huguenin (1971); Backer et al. (1975); Sieber & 
Oster (1975)] — all showing strong fluctuations at a 
frequency of 0.46 cycles per rotation period (here- 
after, c/Pi) and a weaker feature near 0.07 c/Pi. 
Two profile "modes" also have been identified, a 
"B" (for "bright" ) mode in which subpulse "drift" 
is very prominent, and a "Q" (for "quiescent") 
mode wherein the subpulses are weaker and dis- 
organized [Suleymanova & Izvekova (1984); Su- 
leymanova et al. (1998)]. For all these reasons 
the pulsar's subpulse sequences are difficult to ob- 
serve, and our study is based on 430-MHz observa- 
tions benefiting by the high sensitivity of the large 
Arecibo telescope in Puerto Rico. 

A short sequence of "B"-mode pulses is given in 
Figure 1, together with their average profile, and 
the 256-pulse fluctuation spectra in Figure 2 ex- 
hibits both the 0.46 and 0.07 c/Pi features very 
clearly. The former is hardly resolved, whereas 
the latter appears to be just so. A harmonic 
connection between these two features has been 
debated [Taylor & Huguenin (1971); Backer et 
al. (1975); Sieber & Oster (1975)], but never 
clearly determined because of the possible ambi- 
guities involving aliases in an ordinary "longitude- 
resolved" fluctuation spectrum, such as the one 
given in the flgure. We have resolved the mat- 
ter by exploiting the fact that the associated 



modulation is continuously sampled within the fi- 
nite duration of each pulse. Therefore, we first 
Fourier transformed the entire pulse sequence, and 
then, from this imfolded spectrum, established 
that the sets of sidebands associated with the 
features were harmonic. This harmonic-resolved 
spectrum is shown in Figure 3, and will be dis- 
cussed more fully in a subsequent paper (Desh- 
pande & Rankin 1999). It can now be confi- 
dently asserted that the true p/iase-modulation 
frequency is 0.5352±0.0006 c/Pi (though seen at 
0.46 c/Pi as a first-order alias) and that the sec- 
ondary feature at 0.0710±0.0006 c/Pi (a second- 
order alias of 1.07 c/Pi) represents its second har- 
monic. These frequencies vary slightly with time 
and are thus slightly different in our three obser- 
vations, but always fall close to the 1992 October 
values given above. 

With the aliasing resolved, it is also clear that 
the time interval between "drift" bands -usually 
denoted by P3 — is just 1/(0.5352 c/Pi) or some 
1.87 Pi/cycle. Thus, just less than two rotational 
periods are required for subpulse emission to reap- 
pear at the same longitude, both giving the se- 
quence its strong odd-even modulation and estab- 
lishing that the actual, physical subpulse motion 
is from trailing to leading (or right to left) in the 
diagram — that is, negative, or in the direction of 
decreasing rotational longitude. 

The fluctuation spectra of Figs. 2 & 3 also 
show a pair of symmetrical sidebands associated 
with the primary feature. These sidebands fall 
0.027 c/Pi higher and lower than the principal 
fluctuation and represent an amplitude modula- 
tion on the phase modulation. This circumstance, 
viewed along with the remarkable stability of the 
phase fluctuation, demands that these features 
also have an harmonic relationship with the fun- 
damental fluctuation frequency, and so we are not 
surprised to find that (0.535 c/Pi)/(0.027 c/Pi) 
= 20.01±0.08, which is an integer, 20, to well 
within the errors. The entire modulation cycle 
is then just 20 times the 1.87-period one, or some 
37.35±0.02 periods, and this circumstance is beau- 
tifully confirmed in Figure 4, where the full se- 
quence has been folded at this period. 

We have also studied the geometry and polar- 
isation of the pulsar in unprecedented detail. All 
the available mean profiles (seven frequencies be- 
tween 25 and 430 MHz) were used to model the 
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star's conal beam geometry — thus fixing the an- 
gles that the magnetic axis and sightline make 
with the rotation axis, respectively. Four differ- 
ent configurations were explored, corresponding 
to both inner and outer cones [Rankin (1993a,b); 
Mitra & Deshpande (1999)] and to poleward and 
equatorward sightline traverses, respectively. All 
show that at 430 MHz the sightline makes an ex- 
ceedingly tangential traverse through the emis- 
sion cone at close approach still some 102% of 
the half-power radius of the cone. Only a pole- 
ward traverse could be reconciled with both the 
observed polarisation-angle sweep rate of some - 
3°/° (Suleymanova et al. 1998) and the requisite 
magnetic azimuth interval (18°) between adjacent 
subpulses,"'^ which then also fixes the actual, phys- 
ical sense of rotation of the "nearer" (rotation) 
axis as clockwise. We calculate the magnetic az- 
imuth interval between adjacent subpulscs on the 
basis of a) their longitude separation and b) the 
polarisation-angle rotation between them (Rad- 
hakrishnan & Cooke 1969), obtaining values of 
18° independent of (inner or outer) conal radius. 
Clearly, entities at 18° intervals are 20 in number, 
implying a vigesimal beam system circulating in 
magnetic azimuth. 

Pulsar researchers have long tended to view 
"drifting" subpulses as a pattern of "carousel" 
beams [Rudcrman & Sutherland (1975); Sieber & 
Oster (1976)], rotating around the magnetic axis, 
but so far it has not been possible to verify this 
picture either by tracing a given beam through 
its circuit or by establishing the total number of 
beams. Wc sec that our observations and analy- 
sis of pulsar 0943-1-10 now give clear signatures of 
20 beams, rotating counterclock-wise around its 
"emitting" magnetic axis in a total time P3 of 
some 37 periods or 41 seconds. Furthermore, we 
understand the pulsar's emission geometry suffi- 
ciently well that wc can transform the received 
sequence from the observer's frame into a frame 
which is rotating about the magnetic axis of the 
pulsar. This "cartographic" transform then maps 
the intensity as a function of longitude (with re- 
spect to that of the magnetic axis) (j) — (po and 

-'The inncr-conc solution corresponds to a and /3 some 
11.6° and —4.3°, respectively, and all of our calculations 
are based on these values. Other solutions such that 
sin a/ sin/3 Ri —3°/° only change the overall scale of the 
map. 



pulse number k — ko into magnetic-polar colati- 
tude R and azimuth 9 (rotating with period P3). 
Successive pulses then sample the intensity along 
successive chord-like traverses (corresponding to 
different azimuth intervals) through the rotating- 
beam pattern. 

Applying this cartographic transformation to 
the "B"-mode portion of our 1992 observations. 
Figure 5 shows the average configuration of the 
20 subbeams which produce the "drift" sequence. 
The rotation of adjacent beams through our sight- 
line produces the secondary phase modulation, 
and the periodic pattern of varying beam intensi- 
ties gives it its tertiary amplitude modulation. We 
find that the subbeams vary, both in intensity and 
in azimuthal spacing, over 100-second time scales. 
Closely spaced beams tend to be weaker, and some 
beams are observed to bifurcate temporarily, but 
the 20-fold pattern is always maintained strongly 
as a stable configuration. The elongated shape of 
the beams can be understood as a result of trun- 
cation, as most of the emission falls inside of the 
430-MHz sightline traverse. Indeed, a similar map 
at a much lower frequency (111.5 MHz), where we 
expect the sightline to sample emission interior to 
that one above, shows a more nearly circular sub- 
beam shape. 

How sensitive are the conclusions drawn from 
the cartographic transform to any assumptions 
in the analysis? The mapping procedure is 
exquisitely sensitive to the circulation time, and 
only a little less so to the longitude of the magnetic 
axis, to the polarisation-angle sweep rate, and to 
the angles specifying the colatitude of the mag- 
netic axis and sightline "impact" angle. Overall, 
these parameters scale the map — or utterly distort 
it. However, the forward transform used to con- 
struct the maps has a true inverse transform, and 
this inverse cartographic transform can be used to 
"play back" a map in order to produce an artificial 
pulse sequence. Only when this artificial sequence 
correlates in full with the observed sequence do wc 
take the map as correct. Indeed, we have carried 
out such inverse transforms iteratively in a kind 
of "search mode" to refine each of the requisite 
parameters, and a completely consistent picture 
is found in observations from 1974 & 1992 at 430 
MHz and from 1990 at 111 MHz.^ 



^Moving images of these sequences can be 
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The overall structure of 0943+10's radio beam 
provides new insight into the physical nature of 
pulsar emission. Its 20 subbeams rotate almost 
rigidly, maintaining their number and spacing de- 
spite perturbations tending both to bifurcate a 
given beam and to merge adjacent ones — and we 
find an almost identical subbcam configuration at 
the two frequencies. These circumstances suggest 
that the subbeam emission region lies along bun- 
dles of magnetic field lines (or "plasma columns" ) 
having "feet" within a certain annulus on the po- 
lar cap, wherein charges are accelerated at some 
distance above the stellar surface and radiation 
occurs at progressively higher altitudes for lower 
and lower frequencies. We can thus peer at the 
activity at the "feet" of the plasma columns and 
possibly monitor changes viewable at intervals of 
the circulation time. 

The subbeam radiation we observe, estimated 
to be emitted at some 100-300 km height, then re- 
flects processes occurring at lower altitudes along 
these ultra-magnetized plasma columns — indeed, 
some fully down at the surface where the electric 
fields caused by the star's rotation appear. The 
"foot" of each subbeam column moves around the 
magnetic axis, so we cannot appeal to any surface 
features or fixed "hot spots" as their cause. The 
entire rotating-subbeam pattern, averaged over 
time, represents a hollow-cone beam of emission, 
which is emitted along the magnetic-field direc- 
tion at a particular height for a given frequency. 
Signatures of such conical beams are encountered 
in most pulsars, and accurately refiect the angu- 
lar dimension of the magnetic polar cap (Rankin 
1993a). 

Nonetheless, wc have little guidance about 
where the "feet" of these plasma columns will 
fall on the polar cap. Clearly, the radius of the 
active polar cap rpc scales as the square- root 
of the star's radius to that of the light cylinder, 
rpc = er, (27r-r,/cPi)-'^/^ ^ lAbePi^^"^ meters (as- 
suming r* 10 km), and values of e near unity 
or so have been taken to describe the outside edge 
of the emission pattern."^ Therefore, polar-cap 
features corresponding to the subbeam columns 



viewed at sites with the following URL ad- 
dresses: http://www.rri.res.in/~desh and 
http:/ /www. uvm.edu/~jmrankin 
3 e.g., e is (2/3)3/* in Ruderman & Sutherland (1975) and 
1.5 in Rankin (1990), respectively 



would have separations of some 45e meters and 
individual sizes of some 20e m across. These di- 
mensions would increase only a little even if the 
"feet" of the plasma columns happen to be a few 
stellar radii away from the star's surface. Inter- 
preted another way, our estimate of the angular 
width of the subbeams (based on the magnetic 
colatitude spread seen in the map at 111.5 MHz) 
suggests Lorentz factors 7 > 200. We note again 
that our extremely tangential sightline permits us 
to sample only the outermost portion of each sub- 
beam at 430 MHz, and so only this outer portion 
can be mapped down onto the polar cap. 

Ruderman & Sutherland (hereafter, R&S), 
nearly 25 years ago, identified the "drift" -associat- 
ed subbeams with electrical breakdowns (or 
"sparks") in the polar-cap "gap" region, and ar- 
gued that their circulation was due to ExB drift. 
Although the major standpoint of their model 
that even the enormous electric fields that must 
be generated would be inadequate to overcome the 
binding energy of the positive charges (like iron 
ions) at the neutron-star surface — now seems un- 
tenable (Jones 1985, 1986), closely related models 
are still in active discussion [e.g., Usov & Mel- 
rose (1995); Zhang & Qiao (1996); Zhang (1997)]. 
It is noteworthy that no completely independent 
model of subpulse "drift" seems to exist, and our 
maps do suggest a qualitative picture something 
like the model R&S envisioned. Any subpulse 
"drift" model needs to explain first why subpulse- 
scale structures occur and then why they circu- 
late in azimuth. According to Ferraro's Theorem 
(Ferraro & Plumpton 1966) the open-field region 
above the polar cap must rotate with, but per- 
haps more slowly than, the star (so that, for an 
inertial observer, they would both rotate in the 
same direction); thus, in the pulsar's frame the 
regions will rotate oppositely — as is observed for 
0943+10.'' 

Two quantitative features of the R&S model are 
noteworthy. First, their "gap" height is some 50 
m — which is also related in the model to the scale 
over which an active "spark" or plasma column in- 

*That the polar-cap "drift" lags the star's rotational speed 
may, then, be evidence for a "gap" somewhere between 
the emission region and the stellar surface, but it gives no 
indication about whether positive (negative) charges would 
be accelerated outward — or, in R&S's terms, whether the 
star is a pulsar or an antipulsar (Ruderman 1976). 
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hibits the formation of others — and the 45e m we 
infer is fully consistent with this value. Secondly, 
the observed circulation time can also be recon- 
ciled with their model. P3, at polar-cap radius r^, 
for drifting at velocity ExB/S^, can be written 
as P3 = 2-KrsB/<E>Pi, where B and <£■> are 
the magnetic field and the average radial electric 
field in the "gap" region. <E> in turn depends 
on the "gap" potential drop and the "spark" 
distance from the edge of the cap, so that the aver- 
age field is perhaps some <E>= AV/2(rpc — Tg) 
[c./., R&S's eq.(30)]. Assuming that the pulsar's 
magnetic field is correctly computed surface 
value [see de Jager & Nel (1988)'s eq.(l)], then B 
is 2x10^^ Gauss, Pi 1.1 seconds, and thus P3 can 
be as large as 37 Pi only if [rpc — fs) ^ and 
Ay/e^ is well less than 10^^ V. In any case, the 
<E> in the acceleration region is required to be 
^ 3x10^6 Volts/meter, for the circulation to be 
due to ExB drift. 

Our maps suggest that the angular velocity as- 
sociated with the circulation is nearly constant 
across the radial extent of the subbeams, because 
no azimuthal shear is observed within the sub- 
beams. The exact origin and significance of any 
possible dependences of P3 on radius and other pa- 
rameters needs to be understood, and similar esti- 
mates of circulation time in other pulsars (with 
different rotational parameters and geometries) 
should help to define this issue. The finely-tuned 
stability of the pattern observed in this pulsar has 
more to tell. It is unlikely that the remarkably 
periodic and stable arrangement would be possi- 
ble if the required spacing of subbeams were not 
an integral submultiple of the circumference of 
the ring on which they appear to arrange them- 
selves. Pulsar B0943+10's particular characteris- 
tics may therefore result from a critical combina- 
tion of parameters — namely, its rotation period, 
emission geometry, magnetic field strength and 
subbeam spacing, which in turn could be deter- 
mined by the "gap" height. It should not, there- 
fore, be surprising to find that this remarkable 
configuration is indeed an unusual one. 

Nevertheless, averaged over many circulation 
times, 0943-|-10's emission pattern conclusively 
demonstrates the existence of the hollow-conical 
emission beams long attributed to many pulsars 
through less direct means. The subbeam columns 
must be nearly axially symmetric with respect to 



the magnetic axis, because any significant devi- 
ation would require a contrived situation to re- 
produce the periodicities and polarization that we 
observe. Models in which only a part of the polar 
cap is active [e.g., Arons & Scharlemann (1979); 
Arons (1979)] are incompatible with these results. 

The central conclusion to be drawn from our 
subbeam mapping of pulsar B0943-I-10 is that the 
emitting pattern — which is so very stable over 
hundreds of seconds or many circulation cycles — 
is frozen neither to the fields nor to the stellar 
surface. We must then ask where, physically, the 
"mciniory" of the subbeam pattern is carried in 
this system of charges and particle currents de- 
spite the circulation. The detailed structure that 
we have mapped traces the activity at the "feet" 
of the emission columns and represents the first 
direct measurements of some of the parameters — 
such as the locations, size and movement of the 
underlying pattern — reflecting the electrodynam- 
ics not very far from the stellar surface. The over- 
all continuity of the emission, implicit in the stabil- 
ity of the pattern over many circulation periods, 
demands a remarkable "collective steadiness" at 
every stage of the emission process — such as, gen- 
eration and acceleration of particles and coherent 
amplification. 

In summary, pulsar 0943-1-10 exhibits an exqui- 
sitely stable "drifting" -subpulse pattern in its "B" 
mode sequences with a nearly even-odd fluctuation 
frequency of about 0.53 c/Pi. Upon first estab- 
lishing the aliasing-order of this feature, and then 
accurately modeling the star's emission geometry, 
the two "sideband" features were found to result 
from a tertiary amplitude modulation on the sec- 
ondary phase modulation. The remarkable preci- 
sion of the modulation rates argues for a system 
of subbeams circulating around the magnetic axis, 
prompting the use of a novel "cartographic" trans- 
form to map the subbeam structure responsible for 
the pulse-to-pulse fluctuations. The detailed pa- 
rameters of the circulating subbeams allow quan- 
titative assessment of the possible interpretation 
in terms of the system of rotating "sparks" on the 
polar cap as suggested in the Ruderman & Suther- 
land model. The current analysis holds consider- 
able promise for studying the emission properties 
in other pulsars and for assessing physical emission 
theories. 
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Pulse energy Single— pulse sequence 
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Longitude (deg.) 

Fig. 1. — A series of individual pulses from pulsar B0943+10 (centre panel) along with their average (bottom 

panel) and energy (left panel) as a function of pulse number and longitude (360° longitude correspond to one 
stellar rotation). Note the "drifting" subpulses, the alternate-pulse modulation, and the "single" average 
profile. 
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Fig. 2. — Fluctuation-spectral power (centre panel) as a function of longitude and frequency as well as 
the integral spectrum (bottom panel). Note the primary and secondary features at about 0.46 and 0.07 
cycles/period as well as the symmetrical "sidebands" around the former. 
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Fig. 3. — Fluctuation spectrum of a continuously-sampled time sequence reconstructed from the (gated) 
single-pulse sequence used in Fig. 2. The spectrum is presented in a matrix form with each row showing a 
section of the spectrum between n/Pi to (n + 1)/Pi for harmonics n = 0, 1,2, ... of the rotation frequency 
1/Pi. The harmonics at n/Pi are shown separately (left panel) and are essentially the Fourier components 
of the average profile. All other frequency components (up to 100/Pi) are given as a contour plot (centre 
panel), and the column sum of these components, collapsed into a 1/Pi interval (bottom panel), can be 
compared directly with Fig. 2. Note that the principal feature now falls at an (unaliased) frequency of about 
0.535 c/Pi. Note also that whereas the Fourier amplitudes of the principal feature peak at about harmonic 
(row) number 35 (corresponding to a P2 value of some 10.5°), those of the secondary feature peak around 
number 70 — thus demonstrating their harmonicity. 
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Fig. 4. — Pulse sequence folded at the full 37.35-period modulation cycle (centre panel). Twenty emission 
centres can be readily identified, and their relative intensity is shown as a function of longitude (bottom 
panel) and modulation phase in units of pulse number (left panel). For comparison, the dots in the latter 
show the result using a tertiary period that is 2% different than the correct one. 
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Base intensity 



Polar Emission Map 








Radial distance (deg.) 

Fig. 5. — Results of the "cartographic" transformation. The individual-pulse sequence has been mapped 
onto a rotating frame centred on the magnetic axis. Azimuthal average and minimum (unfluctuating "base") 
intensity as a function of magnetic colatitude are given in the bottom and left panels, respectively. Note 
the twenty emission centres and their somewhat elongated shape. A portion of the sightline traverse is also 
indicated at the top of the figure. 
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